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Abstract 

We consider anisotropic cosmological models with an universe of dimension 4 or more, fac- 
torized into n > 2 Ricci-flat spaces, containing an m-component perfect fluid of m non- 
interacting homogeneous minimally coupled scalar fields under special conditions. We de- 
scribe the dynamics of the universe: It has a Kasner-like behaviour near the singularity and 
isotropizes during the expansion to infinity. 

Some of the considered models are integrable, and classical as well as quantum solutions 
are found. Some solutions produce inflation from "nothing". There exist classical asymptot- 
ically anti-de Sitter wormholes, and quantum wormholes with discrete spectrum. 



1. INTRODUCTION 

It is well known that the isotropic cosmological model at present time gives a good de- 
scription of the observables part of the universe. On the other hand, this very fact of our 
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universe's isotropy and homogeneity is puzzling pj. Even in the papers which are devoted 
to the problem of inflation, they start mainly with the metric of the isotropic Friedmann 
universe However, it is possible that at early stages of its evolution the universe exhibits 
an anisotropic behaviour ||. As it was shown in |||, ||], anisotropic cosmological models 
describe the most general approach to the cosmological singularity (the initial singularity at 
some instant to)- Among anisotropic homogeneous models the Kasner solution || represents 
one of the most simple vacuum solutions of the Einstein equations. The Kasner solution is 
defined on a manifold 

M = IR x Mi x M 2 x M 3 , (1) 

where the differentiable manifold Mi (i = 1,2,3) is either IR or S . 

Another very puzzling problem is the fact that the space-time of our universe is 4- 
dimensional. Fashionable theories of unified physical interactions (supergravity or super- 
strings [0, H[ [|) use the Kaluza-Klein idea p0| 1 1 1 of hidden (or extra) dimensions, according 



to which our universe at small (Planckian) scales has a dimension more than four. If the ex- 
tra dimensions are more than a mathematical construct, we should explain what dynamical 
processes lead from a stage with all dimensions developing with the same scale to the actual 
stage of the universe, where we have only four external dimensions and all internal spaces 
have to be compactified and contracted to unobservable scales. 

Exploiting these two remarkable ideas (anisotropy and multidimensionality) , it is natural 
to generalize the manifold (|l|) as follows 



M = IR x Mi x . . . x M, 



n ; 



where Mj (i = 1, . . . ,n) is a dj-dimensional space of constant curvature (or, more generally, 
an Einstein space). If n = 3 and d\ = c?2 = c?3 = 1 or n = 2 and d\ = 2, di = 1 then 
this manifold describes an usual anisotropic homogeneous 4-dimensional space-time. For 
n > 2 and a total dimension D = 1 + J2t=i di > 4 we have an anisotropic multidimensional 
space-time where one of the spaces Mj (say Mi) describes our 3-dimensional external space. 
Multidimensional cosmological models of the type (0) (with arbitrary n) were investigated 



intensively in the recent decade (according to our knowledge, starting from the paper 12] 
investigating the stability of the internal spaces) . 

Quantization of a multidimensional model with a space-time (ph was first performed in 



13]. Beside vacuum models, there were also cosmological models considered which contain 



different types of matter, and exact solutions of the Einstein equations, and of the Wheeler- 



De Witt equations in the quantum case, were obtained (see [14, [L5| and the extended list 
of references there). Exact solutions are of special interest because they can be used for a 
detailed study of evolution of the universe (for example in the approach to the singularity) , 
of the compactification of the internal spaces, and of the behaviour of matter fields. 

In the present paper we consider an anisotropic homogeneous universe of type @/ , where 
all Mi are Ricci-flat. If n = 3 and d\ = di = d% = 1 it describes the usual 4-dimensional 
Bianchi type I model. We investigate this space-time in the presence of m non-interacting 
minimally coupled scalar fields. Scalar fields are now popular in cosmology, because in most 
inflationary models the presence of a scalar field provides homogeneity, isotropy, and almost 
spatial flatness of the universe |2| . It was shown in the paper [ |i~6| ] that for a special form of 
the scalar field potentials these scalar fields are equivalent to a m-component perfect fluid. 
We exploit this equivalence in |T(| to investigate a two-component model (a model with 2 
scalar fields) . Now we shall integrate this model in the presence of 3 scalar fields where one of 
them is equivalent to an ultra-stiff perfect fluid, the second one corresponds to dust, and the 
third one is equivalent to vacuum. The main features of the solutions are the following: If the 
parameters of the model permit the universe to run from the singularity to infinity, then the 
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universe has a Kasner-like behaviour near the singularity, with isotropization when it goes to 
infinity. In the 3-component integrable case, the universe has de Sitter-like behaviour in the 
infinite volume limit. Superficially, it seems this model is not a good candidate for a realistic 
multidimensional cosmology, because of the isotropization of all directions at late times. But 
we shall show that there are particular solutions, which describe a birth of the universe from 
"nothing". The parameters of the model in this case can be chosen in such a way that a 
scale factor of the external space undergoes inflation, while the other scale factors remain 
compactified near Planck length. However this model is really only good, if in addition 
we provide a graceful exit mechanism [jl?]]. For some of the parameters the infinite volume 
limit takes place in the Euclidean region which has asymptotically anti-de Sitter wormhole 
geometry. This is in fact possible because in the infinite volume limit the curvature of the 
spherical spatial sections of the anti-de Sitter geometry decays to zero, whence these spatial 
sections become asymptotically Ricci flat like the spacial sections of our solution. 

Let us also clarify here that, below we refer only to the local properties of (anti-)de Sitter 
space. The relationships between different charts of de Sitter space, with different choices 
of spacial sections, are examined e.g. in IS, 19], while [20| recently provided a classification 
of different multidimensional representations of spaces of constant curvature in arbitrary 
dimension. For a pure d+ 1-dimensional geometry, without additional fields, different choices 
of time, i.e. different slicings into d-dimensional spatial hyper surf aces, should be equivalent 
due to general covariance. Note however that, in [23| and here, the geometry hosts in general 
several additional time-dependent but spatially homogenous matter fields, and the spatial 
homogeneity of any such field is generally not preserved under a change of the slicing. 

Another interesting Euclidean solution represents an instanton which describes tunnelling 
between a Kasner-like universe (a baby universe) and an asymptotically de Sitter universe. 
Sewing a number of these instantons may provide the Coleman mechanism |2l| for the 
vanishing cosmological constant. 

Note that our asymptotically de Sitter solutions are different from the generalized de Sitter 
solutions considered in p^| . There only one factor space (the external one) was Ricci flat, 
and the curvatures of the other (internal) factor spaces were fine tuned with the cosmological 
constant. 

The previous paper [23] has already considered multidimensional cosmological models in 
the presence of a m-component perfect fluid. In the case with one non-Ricci-flat space, say 
Mi, for n = 2 and d\ = 2, d% = 1, this model describes a usual 4-dimensional Kantowski- 
Sachs universe (if M\ has positive constant curvature) or a Bianchi III universe (if M\ has 
negative constant curvature). We also found a 3-component integrable model, where the 
universe has a Kasner-like behaviour near the singularity as in the present paper, but there 
is no isotropization at all. All scale factors corresponding to Ricci-flat factor spaces Mj are 
frozen in the infinite volume limit, but the negative curvature space M\ grows in time. From 
this point of view, the model does not describe usual 4-dimensional space-time, because of 
the missing isotropization, but it may be a good candidate for a multidimensional cosmology, 
if all frozen internal scale factors are near Planck scale. For a positive curvature space Mi, 
the infinite volume limit takes place in the Euclidean region, which there, in contrast to 
the present paper, has wormhole geometry only w.r.t. the space Mi, and the wormhole is 
asymptotically flat. 

In the present paper we consider homogeneous minimally coupled scalar fields as a matter 
source. Usually, only real scalar fields are taken. Here we admit also purely imaginary scalar 
fields. Such scalar fields imply a negative sign at the kinetic term in the Lagrangian. Such 
scalar fields may arise after conformal transformation of real scalar fields with arbitrary 



coupling to gravity [24, 25, p7y . They appear also in the Brans-Dicke theories after the 
dimensional reduction from higher dimensional theories H |28|, [2j| . Also the C-field of Hoyle 
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and Narlikar has a negative sign in front of the kinetic term pC| ]. The authors of [31, 32] 
emphasize the need for scalar fields with negative kinetic terms in multidimensional theories 
in order to fit the observable data (see also a discussion of this topic in |33|]). As we will 
show here, in the particular case of constant if, the imaginary scalar field is equivalent to a 
negative cosmological constant which results in an anti-de Sitter universe. In what follows 
we do not exclude the possible existence of imaginary scalar fields, whence in our paper we 
consider real as well as imaginary scalar fields. 

This paper is organized as follows. In section 2 we describe our model and get an effec- 
tive perfect fluid Lagrangian, exploiting the equivalence between an m-component perfect 
fluid and m non-interacting scalar fields with a special class of potentials. In section 3, we 
investigate the general dynamics of the universe and its asymptotic behaviour. In section 4, 
classical solutions for the integrable 3-component models are obtained. Classical wormhole 
solutions are obtained in section 5 where it is also shown that they are asymptotically anti-de 
Sitter wormholes. Section 6 is devoted to the reconstruction of the scalar field potentials. 
Solutions to the quantized models are presented in section 7. In section 8 we summarize our 
results. 



2. THE MODEL 

We consider a cosmological model on a multidimensional manifold (||) with a metric 



g = g MN dx M ® dx N = -e 2 ^dr ® dr + £ e 2 ^g®, 



(3) 



where, for i = 1, . . . , n, = g m \ ni dx mx (g> dx n \ m,, = 1, . . . ,<£*, is the metric form of the 
Ricci-flat factor space Mj dimension di. 



R 



mini 



0, i = 1, . . . n. 



The action of the model is taken in the form 



S = J d D xJ\g~\R[g] +S ip + S l 



(4) 



(5) 



where <S GH is the standard Gibbons-Hawking boundary term, k 2 is the gravitational coupling 
constant in dimension D = Y%=i di+1, and = J2T=i is the action of m non- interacting 
minimally coupled homogeneous scalar fields 



SjO = - / d D xJ\g] \g MN d M V> {a) d N ip^ + U^ a \ip^) 



For the metric (|3|) the action (||) reads 



5 = fj, J drL s , 



with the Lagrangian 

L s = -e 
2 



-7+70 



>h2 



2 e 7+70 ^ U {a) (ip (a) ] 



(6) 



(7) 



(8) 



a=l 



a=l 
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Here 70 = Y%=1 dif3 l and \x = rii=i Vil n where Vi is the volume of the finite Ricci-flat spaces 
(Mi,g^). The components of the minisuperspace metric read 

Gij = di5ij - didj. (9) 

As in [16 1 we subject the scalar fields to the perfect fluid energy- momentum constraints 

p(«) = ( Q ( fl ) - 1) p ( a \ (10) 

with constants = 1, . . . , m, and the energy densities 

p (a) = _ T (a)0 = 1 -2 7( ^( a ))2 + ^(a)^(«)) ( n ) 

and momenta 



p(a) = T (a)JJ = 1-27 ^(a))2_ ^(0)^(0)^ M = 1, ...,£) — 1 , (12) 



according to the Lagrangian (|8|). In 16] it was proved that, for cosmological models with 
a metric (|3|), the presence of m non-interacting scalar fields satisfying the relations (pX|) is 
equivalent to the presence of an m-component perfect fluid with a Lagrangian 

1 m 

L p = - e -^°G ij $ i $ j - K 2 e 7+70 2 /) (a) , (13) 

2 a=l 

and energy densities of the form 

p (a) =yl (a) r o(«' i a=l,...,m, (14) 
with constants A( a ) and a spatial volume scale 

n 

V = e 70 = J| af (15) 
i=l 

defined by the scale factors Oj = e' 3 ', i = 1, . . . ,n. Note, that the total spatial volume is 
Vtot = fi'V. The energy density p( a ) is then connected with the pressure via fllOD, and 
equations ( pT| ) and ( |l2| ) imply a*-")//' 1 ) = e~ 27 (i^ a )) 2 . So, for real scalar fields and positive 
a.( a \ the energy density of the perfect fluid is positive. But, keeping in mind the possibility 
of imaginary scalar fields (see Introduction), for the general model we shall also consider the 
case p( a ) < 0. Then, the constants may have any sign. 

Assuming the speed of sound in each component of the perfect fluid to be less than the 
speed of light, 

-\p {a) \ < P (a) < \p {a) \, a = l,...,m. (16) 



With ( |T0| ) this implies the inequalities 

< a (a) < 2, a = l,...,m. (17) 

Note that, with p = £™ =1 p(°) and P = £™ =1 P {a \ the energy dominance condition requires 
only — \p\ < P < |p|, rather than (|l^). In this paper however, although it might be 
possible to generalize results for arbitrary a^ a \ for simplicity we keep the assumption ( |T6| ) 
in order to make use of the inequalities ([17]) . 
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Exploiting the mentioned equivalence between scalar fields and perfect fluid, we investi- 
gate the dynamics of the universe via the Euler-Lagrange equations of (|l3|), and reconstruct 
the scalar field potentials f7(°)(y>W) satisfying the perfect fluid constraint (|l0|) . 



3. GENERAL DYNAMICS OF THE UNIVERSE 

In the harmonic time gauge 7 = 70 = J2?=i di(3 l (see e.g. [13, |25|), the Lagrangian (0) with 
energy densities ( |l4] ) just reads 

-. m 

L^-GaPP -«V^/ a ». (18) 

1 0=1 

Then the corresponding scalar (zero energy) constraint can be imposed as 

-. m 

-Gijftfr + kV 70 ^ p (a) = 0. (19) 
2 ' 0=1 

The minisuperspace metric may be diagonalized (see also |l3|]) to 

n-l 

G = r/ fc/ dz fc (8) dz' = -dz° (8) dz° + ^ <te' (8) (20) 

where 



i=l 



z° 



z l 



% = 1, . . . , n — 1, and 



? _1 E^' =<?" 1 ln^ (21) 
i=l 

n 

[di/Ei^+i] 1 / 2 £ dj , (22) 



g:=[( J D-l)/( J D-2)] 1 / 2 , (23) 

i=k 

With the aid of these transformations the Lagrangian (18) and the scalar constraint ( |l9| ) can 
be rewritten as 

-, m 

L p = - Vkl z k z l -k 2 Y, exp(^V), (24) 
1 0=1 

-77 W i fc i' + k 2 E ^ (a) exp(fe (a) gz°) = (25) 



a=l 



respectively. Here, k^ := 2 — a^ a ^ (a = 1, . . . , m), whence the inequalities ([[?]) for hold 
also for /c( a \ 

< k^ < 2. (26) 
The equations of motion for z l , i = l,...,n — 1, simply read 

t = 0, (27) 
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and readily yield 



z l = pW + q\ 



(28) 



where r is the harmonic time, p l and q l are constants. Clearly the geometry is real if p l and 
q % are real. The dynamics of z° is then given by the scalar constraint (25), which may now 
be written as 

1 



i0\2 



+ e + k 2 A(a) eMk (a) qz°) = 0, 



a=l 



for a real geometry with 



1 



n-l 



i=l 



The coordinate transformations ( pip and (22) can be written as 

n 

z k =Y J t k i(i i , k = 0,...,n-l, 



(29) 



(30) 



(31) 



i=l 



whence the inverse is given by 



n-l 



/3 l = E* 



% = 1, . . . , n. 



k=0 



For % = 1, . . . , n, with i°j = di/q and t l Q = [q(D — 2)] 1 we obtain the scale factors 

a i = A i V 1 ^ D - 1 h aiT , 



where 



n-l 



n—l 



i=l 
n-l 

i=l 1=1 



l=i l=i 

The parameters a 1 satisfy the relations 

n 

E = 0) 

i=l 

and, analogously the parameters 7* fulfil 

n 

E^y = 0, 

i=l 

n n—l 

E^(V) 2 = E(^) 2 - 

i=l 1=1 
From the definition ( |34| ) and the relation (37) it follows that 

f[Af = 1. 



(32) 

(33) 
(34) 

(35) 
(36) 



(37) 
(38) 



(39) 



Note that using the constraints (|35|) and (|39|) the equation (|3) yields again (15), i.e. 



11^=1 a f = V • Recall that r in ( |33|) is the harmonic time. The synchronous and harmonic 
times are related by 

t = ±Je" /o dT + t = ±JvdT + t . (40) 



The expression ( J33| ) shows that the general model does not belong to a class with static 
internal spaces (see e.g. |34|]), but just for e = ( i.e. a 1 = 0, i = l,...,n), there is a 
solution 

a i = A i V 1 ^ D ~ 1 \ i = l,...,n, (41) 

which is isotropic. 

In order to find the dynamical behaviour of the universe we should now solve the con- 
straint (|29|), i.e. the mechanical energy conservation equation 



e = T + U (42) 

with kinetic energy T := ^(i ) 2 and potential U := — k 2 J2™=i exp(k ( - a " 1 qz°). Depending 
on the parameters and their signs, the potential U may exhibit a rich structure with 
several extrema, and a classical Lorentzian trajectory is bound by possible turning points at 
e = U. Since the general dynamics is very complex, we investigate the asymptotic behaviour 
of our model universe in the limit of large spatial geometries V — ► oo and near the singularity 
V —* 0. Without restriction we suppose now 

< fcW < . . . < fc (m) < 2. (43) 

1. Limit V — > oo: In the limit V — > oo (i.e. zo — ► oo) the term —k 2 A^ exp(k^ m \z°) domi- 
nates the potential U, whence, for k^ ^ 0, there are two cases to be distinguished: 

i) ^4.( m ) > 0: Here, for V — ► oo, the term e may be neglected. So, the constraint equation 



( f42|) has the asymptotic solution 

e«*> = V w (2K 2 AM)-Vk Cm) (g| T |)-2/* (m) , (44) 

with 2g := k^q, where (without restriction) we have chosen initial conditions such that 
V — > oo at r — > 0, when according to ( |33"1) the system is subject to an isotropization, 

a^y 1 /^" 1 ), V^oo, i = l,...,n. (45) 

In this limit the harmonic and synchronous times are connected by 

| r | „ | t |fcW/(fe(™)-2) ; fc (m)^ 2> (46) 



t ~ exp 



(-V2^lM 9 |t|) , fc (m) = 2 . (47) 



So the synchronous time evolution of the spatial volume is (asymptotically for t — ► oo) given 
by 

y ~ \t\ 2 ^ (m) , k^^2, (48) 
F ~ exp(y/2K 2 A( m )q\t\), k (m) = 2 , (49) 
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with scale factors (according to isotropization) 

ui2/a("0(U-l) ^ k (m) ^ 2 



04 ~ 



Qi~exp( — - |t|) , fcW=2. 



(50) 
(51) 



Taking a usual anisotropic space-time model {D = 4, n = 3, d\ = c?2 = c?3 = 1) then for large 
(synchronous) times the formulas ( |50| ) and ( |5l| ) yield scale factors a% ~ |i| 2//3 for /c^™- 1 = 1 

(dust) and dj ~ exp (\J fc 2 ^4( m )/3|i|) for k^ = 2 (vacuum). Asymptotically, power-law 

inflation (with power p > 1) takes place for < a( m ) < 2/(D - 1), and = 2/(D - 1) 
yields a generalized Milne universe. 

ii) j4( m ) < 0: Here, the Lorentzian region has a boundary at the turning point V mSuX of 
the volume scale, which in the large energy limit e — > oo is asymptotically given as 



* rn 



K 2 |A( m )| 



(52) 



The region with V > V ma _ x is the Euclidean sector. For V >> V max , we obtain the asymp- 
totically isotropic solution 



2K 2 |A( m )|g|r| 



2/fc( m '(D-l) 



(53) 



In the Euclidean region, we obtain a classical wormhole w.r.t. each factor space. With 
constants of integration (in (]28|)) pi = (i.e. a* = 0), i = 1, . . . ,n, the wormhole takes its 
most simple and symmetric form. Then the throats are given by 



'(th)i 



m) I 



(54) 



In the case k^ = 2 we obtain asymptotically (for t — * 00) anti-de Sitter wormholes with 
synchronous time scale factors 



l2K 2 \A( m )\q 
ai ~ exp ( v |t| 



,n . 



(55) 



2. Limit V — » 0: For fcW 7^ 0, in the small volume limit V — > 0, i.e. z° — > —00, the potential 
vanishes £/ — ► 0. So, for e > 0, we obtain (asymptotically for t — > 0) a (multidimensional) 
Kasner universe [Ell, |36|1 , with scale factors 



a,; ~ i 



1, . . . ,n . 



with parameters a 1 satisfying 



(56) 



(57) 



i=l 



i=l 



If fc(!) = 0, then U -> -k 2 ^ 1 ) for z° -> -00. Here, for £ := e + « 2 AW > 0, we obtain 
(asymptotically for t — > 0) a generalized Kasner universe [35|, i.e. scale factors fl5"q ) with 
parameters a* satisfying 



J^a* = 1 , J^dj 
i=l i=i 



a') 2 = 1 - a 2 , 



(58) 
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with the parameter a — ► for — ► 0. 

In the exceptional case = e + k 2 ^ 1 ) = the term of the matter component a 
dominates the constraint (|29|) , whence we obtain (compare also fl6TQ 



v ~ t 2 / Q(2) , 

, t 2/[(D-l)aW ]exp f i /(a (2) )t -^y 



a 



,n , 



(59) 
(60) 



where f(x) := 



(2-£) 



(2-a;) 2 ( 3 2 K 2 A( 2 ) 



In another exceptional case where e = (i.e. 
a" = 0, i = 1, . . . , n) the universe is isotropic everywhere, i.e. a% ~ 

yl/D-l i i = i }mmmt n. If, 
for example, k^ = (and = 0) we obtain from (|59|) or (BO) 



a,- ~ t 



2/[(D-l)a( 2 )] 



(61) 



4. INTEGRABLE 3-COMPONENT MODEL. CLASSICAL SOLUTIONS 

In this chapter, we consider the integrable case of a three-component perfect fluid (m = 3) 
where one of them (a = 1) is ultra-stiff matter (fcW = 0,a^ = 2), the second one (a = 2) is 
dust (fc( 2 ) = l,a( 2 -* = 1), and the third one (a = 3) is vacuum (k^ = 2,a^ = 0). The case 
jfe(i) = 0, fc( 3 ) = 2fc( 2 ) with < < 2 is also included if one substitutes q by q = k^ 2 'q. 
The constraint equation (p9|) reads in this case 



2 

and can be rewritten like 



where 

and the potential U(z°) is 



I (>) 2 + £ + K 2 A {X) + K 2 A (2) e ^ + ^2^(3)^0 = ( g 2) 

i(i°) 2 + ^(z°), (63) 
£ := e + k 2 ,4« (64) 



E 

2 



I7(je ) := -Be qz ° - Ce 2qz ° (65) 

with the definitions B := k 2 A^> and C := k 2 ^ 3 ). 

As mentioned in the introduction, for a complete description of the model the parame- 
ters E, B, and C are considered to take positive and negative values. Then, we have four 
qualitatively different shapes of the potential ( pl|) (see Fig. 1 and Fig. 2). For each of them, 
we shall solve the constraint equation separately. Eq. (|63] ) integrates to 



±V2q(r-T ), (66) 



WE + BV + CV 2 
where r is the harmonic time coordinate, and tq is a constant of integration 
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i) B > 0, C > (see Fig. 1): The solutions of equations (36) are 



V 



V 



B 3l 

2 



f 2 _ 2C 
J 

4Ef 



(B - /) 2 -4EC" 
2|£7| 



-B = 0, 

£ > 0, 
E < 0, 



(67) 

(68) 
(69) 



where A := AEC - B 2 (|A| = B 2 + 4|£|C for £ < 0) and 



f = T-T , E = 0, 



^2C 



qB 



< It — Tq I < +00 



/ = exp (V2Eq{r - r )) , £ > 0, In (b + 2\ZIc) < In / < +oo , 
/ = sin (\j2\E\q(r — tq) \ , E < 0, — arcsin ( —7== ) < arcsin / < 



IAI 



The synchronous and harmonic time coordinate are related via 



(70) 
(71) 

(72) 



r - r = coth I \j ^-qt | , £ = 0, 



/C 



exp (V2Eq{T - r ) ) = 5 + V4£ , C coth ( J— gt ] , £ > 0, 



tan J|B|/2g(T-T ; 



IAI 



B 



^cothfJ^i-J 



£ < 0. 



(73) 
(74) 
(75) 



Using these relations, we obtain the expressions for the volume of the universe in synchronous 
time: 



B 

V = — sinh 



C 



V 



c 



B + y/AEC coth I \l-qt 



sinh 



2 

IC 



</l 1 . E = 0, |i| < 00, (76) 
gt I , -E > 0, < t < +00, (77) 



V 



2|£|(l + tan 2 (y/2)) 



5(1 + tan 2 (y/2)) + 2y^\A\ tan(y/2) 



, £<0, 



(78) 



where tan(y/2) = tan K/|.E|/2g(T — To) J is given by ([7g). Expression ([Tq ) can be written in 
a more convenient way if the parameter To is chosen such that equation (|69|) is symmetric 
with respect to the turning point Vq = (^—B + y/\A\J /2C, namely 



V 



2\E\ 



t < 



B + v/|A[cos (v/2|B|gT) ' V 2 \ e \Q . 



IT IB 

h arcsin — === 

2 \ ,AN 



(79) 
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In this case, 



tan J\E\/2qr 



v/lAT-B 



tanh 



-qt 



(80) 



and for the volume results 
1 



V 



2C 



A \-B + 



|A| + B ) tanh 2 | \/— 



cosh 



|i|<oo. (81) 



The region V < Vb is the Euclidean sector and we obtain the instanton by analytic continu- 
ation t — > —it in formula (j8l|): 



V 



2C 



B 



|A| ! .0 ] Ian 2 | \l^qt 



cos 



-qt 



(82) 



with \t\ < ^j- arctan . / V^i — ■ 

On the quantum level, this instanton is responsible for the birth of the universe from 
"nothing" . 

ii) B < 0, C > (see Fig. 2): In this case, the maximal value of the potential U(z°) is U m = 
B 2 /AC at = - In (\B\/2C) and for < E < V m we have two turning points, namely: 



V,- 



(1,2) 



B\ ± J\A\ /2C 



(83) 



(i) 



where |A| = B — AEC '. Classical motion takes place either for < V < Vq or for 

y (2) < v < +oo. 

If E 1 < 0, we have one turning point only, namely 

V = \B\/C, E = 0, 



(84) 
(85) 



B\ + J\A\) /2C, E<0. 



where |A| = B 2 + 4|i?|C and classical motion takes place for V > Vq. 
The solutions of the equation (|66|) read 



V = 

V 
V 
V 



1 



\B\ Q 2 ( 2C 
2 \B 2 q' 2 


-P) 


AEf 




(\B\+f) 2 ~ 


AEC 


AEf 




AEC - (\B\ 


-f) 2 


AEf 





E = 0, 



0<E<U m , 0<V<V c 



(i) 



o > 



(86) 
(87) 



< E < U m , V^ 2) < V < +oo, 



(\B\ +f) 2 -AEC 



V 



2\E\ 



W\f~\B\ 



E>U m , 
E < 0, 



(90) 
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where 



/ = r - T , E = 0, 

\t-t \ < V2C/q\B\ 



f = exp (V2Eq(T - r 
/ = exp (V2Eq(T - r 
/ = exp (V2Eq(T - r 



0<E<U m , V <Vr 



(i) 



' 



hW|A| < In/ < +oo, 
< E < U m , V > Vr 



(2) 
' 



/ = sin J2\E\q(r - t ] 



ln(\B\ - V^ECj < In / < \nyf\A\ 

E>U m , 

ln(-|B| + V4£?c) <ln/<+oo, 
£ < 0, 

\B\ 



(91) 
(92) 

(93) 

(94) 

(95) 



7T 



arcsm ■ 



< arcsin f < —. 
|A| - 2 



The harmonic and synchronous time coordinates are related via 



tan J\E\/2qr 



2C 

t = | — ;— tanh 
\B\q 

f = V4EC cot 
f = -V4EC tanh 

f = v / 4ECcoth 

^/4\E\C 
J\A\ + \B\ 



-qt 



(96) 



2 
2 



E = 0, 

qt]-\B\, 0<E<B 2 /4C, V < V Q {1) , (97) 
qt]+\B\, 0<E<B 2 /AC, V > Vq 2 \ (98) 



2 

tanh 



'// I " 1^1 

C 
-qt 



E > B 2 /4C, 



E < 0, 



(99) 



(100) 



where in (|9q ) and ( pLOOj) the constant ro is chosen such that the expressions are symmetric 
with respect to the turning point at the minimum. Then, the volume of the universe is 



V 



V 



\B\ 



IC 



i^cosh z ( \l-qt | , E = 0, \t\ < +oo, 



1 

C 



IC 



V4EC coth \ \l-qt ] - \B\ 



< E < B 2 /AC, V <V^\ 0<t< 



sinlr j \l^qt \ . 



(101) 
(102) 



'2Cq 



-arcoth 



V = — 



1 

2C 



\B\ + \/|A| - j \D\ - \/|A| ) lank- j 



|£| + V]A[ 



cosh I ^/ — 



(103) 
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0<E<B 2 /AC, V>V {2 \ \t\<+oo 



V 



1 

c 



\/4ECcoth \\l-qt 




\B\ 



sinh 




E > B 2 /AC, < t < +00 



V = — 



1 

2C 



A| + \B\ + [ J\A\ - \B\ tanh 2 \ h-qt 




(104) 



(105) 



E < 0, |i| < +00. 



Eqs. (101), ( |103j ), and ( |105| ) are written in a symmetric way with respect to the turning 
point at t = 0. The instanton solutions can be obtained by analytic continuation of these 
symmetric expressions and result in: 



V 



V 



V 



cos 




\B\ 
~~C 

1 

2C 



0<E<B 2 /AC (|A| = 5 2 -4EC), \t\ < vr/V2Cg, 



qt\ , E = 0, \t\< n/V2Cq, 



A I tan 2 



B\ + ,/IAI + \B 




1 

2C 



| J B| + ,/|A|+ \B 



A I tan 2 




cos 




(106) 
(107) 

(108) 



E<0 (|A| = B 2 +4\E\C), \t\< 



>2Cq 



■ arctan , 



\VW\-\b\ 



( 2) 

In equation ( |107| ), the instanton is symmetric with respect to the turning point Vq . For 
the same instanton but now symmetric with respect to the turning point at Vq 1 ^ , we have 



V 



1 

2C 



\B\ 



A| + 151 + \/|A| tan 




cos 




(109) 



< E < B 2 /4C, \t\ < n/V2Cq. 



All the instantons ( |106| ) to (|l09j) are responsible on the quantum level for the birth of the uni- 
verse from "nothing". The instanton ( |107[ ), ( |109| ) is of special interest. Its qualitative shape 



is seen in Fig. 3 where V m i n = Vg and V max = Vq . The instanton describes tunnelling 
between a multidimensional Kasner-like universe (a baby universe) and a multidimensional 
de Sitter universe because, as was mentioned in chapter 3 and as we shall demonstrate more 
precisely latter, the limit V — > corresponds to a Kasner-like universe and in the limit 
V — > 00 we obtain an (isotropic) de Sitter universe (in [37] to [40] analogous types of an 
instanton describing tunnelling between a Friedmann universe and a de Sitter universe were 
obtained for a different model) . The instanton may also represent the birth (demise) of a de 
Sitter universe (see ( |107| )) and a baby universe (see ( |109| )) from (into) "nothing". As was 
demonstrated in p7|, |38|], the instanton may be extended beyond V = V m i n ,max gluing to- 
gether a number of Euclidean manifolds (see Fig. 4). Such gluing may provide the Coleman 
mechanism [21] that establishes the vanishing of the cosmological constant. 
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The case with E = B 2 /AC is degenerated. Here, two turning points coincide with each 



other: V c 



(l) 



o 



T/( 2 ) - V 



\B\/2C. In this case, we obtain with the synchronous time 



coordinate t for the two volumes the expressions 

LB I 



v -m 



1 — exp 



'2E 



qt 



< t < +oo, 



(110) 



which describes an infinitely long lasting rolling down from the unstable equilibrium position 
Vq to the singularity V = and 



V 



2E 

W\ 



1 + exp 



\B\ 



'2E 



qt 



\t\ < +oo 



(111) 



describing the infinitely long lasting rolling down with V — > oo. 

To obtain solutions in the two remaining cases B > 0, C < (see Fig. 2) and 
B < 0, C < (see Fig. 1), it is not necessary to solve equation ( |66| ) again. We can 
instead take the solutions found already in subsections i) and ii). It is clear that the Eu- 
clidean solutions obtained there are Lorentzian ones here and vice versa Lorentzian solutions 
of subsections i) and ii) are Euclidean ones here. What we have to do is the evident substi- 
tutions B — > \B\, C — > |C|, and E — > —E where it is necessary. For example: 



iii) B > 0, C < 0: From fllQ6p , ( gOp , and ( jjjg ), we obtain respectively 
V 



— cos 2 ( V|C|/2^), £ = 0, \t\<n/ qx /2\C\, 



1 



A| + £ 



V 



2\C\ 
E > 0, |*| < 
1 



|A| - B) tan 2 U/|C|/2gt) 
|A|+^/ 



cos 



arctan 



A\ — B 



\C\/2qt 
li/2 



(112) 
(113) 



A| + LB + */|A| tan^ J\C\/2qt 



2|C| 

-B 2 /4|C| < £ < 0, |t| < n/J2\C\q . 



cos' [J\C\/2qt , (114) 



Solution ( |113| ) is symmetric with respect to the classical turning point. To investigate the 
limit |t| — ► it is better to give another representation of the same solution, namely 



V 



1 

W\ 



B + ^4E\C\cot ( \/\C\/2qt 
E>0, < * < 



BW*[J\C\/2qt 



(115) 



2 /A"-fi 
arctan ■ 



y%C\q ^ V4E\C\ 
In this case, the harmonic time coordinate and the synchronous one are related via 



exp (V2Eq(T -t ))=B + JAE\C\ cot J\C\/2qt 



(116) 

Equation ( |114| ) is symmetric with respect to the turning point Vg 1 ^ = (b — y^AM /2\C\. Its 
analytic continuation gives a "parent instanton" (see ( |102|) ) with 

V = — \V4EC coth U\C\/2qt\ - B sinh 2 (j\C\/2qt\ , (117) 



\C\ 



-B Z /4\C\ <E<0, 0<t< 



-arcoth- 



V4EC ' 
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which is responsible for the birth of a baby universe from "nothing" . The Lorentzian solution 
(114) symmetrically written with respect to the turning point Vq 2 ^ = (^B + y^Ajj /2|C| reads 
(see (p7D ) 

V 



1 



B + J\A\ + B 



2|C| 

-B 2 /4|C| < £ < 0, \t\ < n/j2\C\q 



A| tan 2 J\C\/2qt 



C os 2 (J\C\/2qt), (118) 



we obtain 



iv) B < 0, C < 0: Here, a Lorentzian region exists for E > only. Prom equation (|66|), 

V = ^[^ICIcot^V^^-lSljsin 2 ^^^^, (119) 

2 



E > 0, < i < 



-arccot- 



sin^ I V|C|/2^ 
|A| + l-BI 



and the equation relating the harmonic time coordinate and the synchronous one reads 



cxp 



(V2Eq( 



r - T , 



4E|C|cot J\C\/2qt) - \B\. 



(120) 



These equations are useful for the investigation of the small time limit |t| — > 0. 

To obtain an instanton solution (wormhole) it is necessary to rewrite equation ( |119| ) 

symmetrically with respect to the classical turning point Vq = (-\B\ + VlA|) /2\C\. We 

can reformulate equation ( |119| ) or use directly equation (82). The result is 



V 



1 



2\C\ 
E>0 



\B\ 



A| + \B\ tan 2 J\C\/2qt) 



cos z [J\C\/2qt), (121) 



t < 



VW\q 



arctan , 



IAI - \B\ 



We shall investigate now the small time limit \t\ — ► for Lorentzian solutions. As we 
shall see, it corresponds to the vanishing volume limit V — > and takes place for E > if 
B, C > or B > 0, C < and for E > if B, C < or B < 0, C > (see Fig. 1 and Fig. 
2). First, we consider the case of positive energies E > 0. As follows from @, @), (PD , 
and fll20D 



exp ( v2Eqr ) ~ t 







(122) 



and with the help of equations (^), (^), and ( p6[ ) we obtain for the scale factors in this 
limit the expressions 

ai^Ait & \ i->0, (123) 

where 



1 



1 



o 



D- 1 

and the parameters satisfy the conditions 

1=1 



7 2^g 



-Q 



(124) 



(125) 
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x>(«*) 2 = 1 



i=X 



1 k 2 AW 
q 2 e + k 2 AW 



1 for 







(126) 



in accordance with the Eqs. (pq) to (pq). For the volume of the universe, we obtain in this 
limit 

V ~ t, t^O. (127) 

Thus, for positive energy, E > 0, and small synchronous times the universe behaves like the 
Kasner universe. 

Now, we consider the exceptional case E = 0. It follows from (|73|) that 



2 1 



and for the volume, we obtain from (|7 



q 2 Bt 



V~t 2 , t 



t -> 



0. 



(128) 



(129) 



With the help of Eq. (|33| ) , we conclude that the approximation of the scale factors is given 
by 

/2a; r 



a^i^ 13 - 1 ) exp 



W 2 £t/ ' 



t -> 



(130) 



in accordance with expression ph for a( 2 ) = 1. 

Thus, the scale factors behave either anisotropically and exponentially like 



exp 



2a, ; 1 







.q 2 BtJ' 

if aj 7^ (e > 0, yl'- 1 ) < 0) or they have power law behaviour like 



(131) 



(132) 



t 2 , t 



if a, = (e = 0, = 0) (see Eq. ( |6l|) for a^ 2 ^ = 1). In the latter case, the free minimally 
coupled scalar field is absent (A^ = 0). 

Similar investigations can be done for the equation (112) shifted in time such that V ~ 
0. 

If E < 0, the universe has in the Lorentzian region a classical turning point at the minimal 
volume V m i n and reaches never V = (see Fig. 1 and Fig. 2). 

Now, let us consider the infinite volume limit V — ► which, as we shall see, coincides 
with the limit t — ► +oo. As follows from Fig. 1 and Fig. 2, this is possible if B, C > or 
B < 0, C > 0. With Eqs. (|73|) to (|7|) and (!§, @ to fllOOl) one can demonstrate that 
r becomes asymptotically constant for t — > +oo and the constant can be put equal to zero 
(with a proper choice of the integration constant tq). From (^) follows that isotropization 
takes place in this limit, namely 



yl/D-X 



+oo 



and from Eqs. (|7§) to © and (|i0l| ), ( gOp to ( |105D we get 

V w exp (\/2Cgi) , t +cxd 



(133) 



(134) 



in accordance with (4£). 
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Thus, if C > we obtain in the limit t — ► +00 an (isotropic) de Sitter universe. If C < 0, 
the universe has a classical turning point at maximal volume V max and the volume can not 
go to infinity. 

Let us come back once more to the case C > describing a universe arising from "noth- 
ing". The volume is given by (|8l]), (|105|), and ( 103| ) and the harmonic time coordinate and 



the synchronous one are related via fl80D, ( |10C ), and fl9q) , respectively. We shall restrict 
ourselves to the case E < for simplicity. In this case, we get the asymptotic expression 

arctan -= \ — = A, (135) 



if t S> (h is sufficient to satisfy \J~^qt > 2). Then, as follows from equation (p3|), 

the scale factors are given by 

Oi « ^ exp(aM)F 1/(D - 1) . (136) 



In [17 1 was shown that for 4 $ y^rqt <C -D — 1 the parameters of the model can be chosen 



2' 

such that, due to the exponential prefactor in ( |136|) , some of the factor spaces (with a % > 
0) undergo inflation after birth from "nothing" while other factor spaces (with a' 1 < 0) 
remain compactified near the Planck length L P1 ~ 10~ 33 cm. The (graceful exit) mechanism 
responsible for the transition from the inflationary stage to the Kasner-like stage, in which the 
scale factors of the external spaces Mj exhibit power-law behaviour while the internal spaces 
remain frozen in near the Planck scale, deserves still more detailed investigations, similar 
those of J!?]]. (There the complementary case of multidimensional cosmological models with 
cosmological constant was considered.) 



5. CLASSICAL WORMHOLES 



In this chapter we consider in more detail a special type of instantons, called wormholes. 
These usually are Riemannian metrics, consisting of two large regions joined by a narrow 
throat (handle). Obviously, they appear if the classical Lorentzian solutions of the model 
have turning points at some maximum, Namely, according to Fig. 1 and 2, for models with 
C < (the parameter B may be positive as well as negative). Let us show this explicitly. We 
consider instantons which can be obtained by analytic continuation t — > — it of the Lorentzian 
solutions ( |TT2|) , (|lia| ), O and (|l2l|) respectively. 



V 
V 

V 
V 



cosh z (J\C\/2qt , E = 0, \t\ < 00, 



B 

W\ 

1 



2\C\ 

E > 0, \t\ < 00, 



A| +B+ [J\A\ -B tanh 2 [J\C\/2qt) cosh 2 (J\C\/2qt , (138) 



1 



2\C\ 

-B 2 /A\C\ < E < 0, |t| < 00 



B + J\A\- [B-J\A\ )t&nh 2 (J\C\/2qt) cosh 2 (J\C\/2qt ), (139) 



1 



2\C\ 

E > 0, \t\ < 00. 



A| - \B\ + J\A\ + \B\ tanh 



\C\/2qt 



cosh 2 [J\C\/2qt\, (140) 



(137) 
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As mentioned before, these equations correspond (with evident substitutions) to the Lorentzian 
equations ( |101[ ), ( |105| ), ( |103j ) and (81), respectively. The harmonic and synchronous times 
are related respectively by 



Bq 



tanh J\C\/2qt 



tan \ E/2qr 



tanh J\E\/2qr 



s/\A\+B 
v/IST-B 



tanh \\C\/2qt 



ViEC 



tanh J\C\/2qt) 



E = 0, 

E>0, 
B 2 



tan (j/E/2<lS) = tanh (j/\C\/2qt) 



4|C| 
E > 



< E < 0, 



(141) 
(142) 
(143) 
(144) 



(See d§, ( |i00|) , (H) and (|0]) respectively. (||) looks like M), if we choose the constant of 
integration ro such that /| r =o = ^o^' whence / = \/\A\ exp ^\/2-Egr^ , and use the relation 

.|b|-vTaT 



/ = \B\-V4EC tanh 



y / C/2qt + artanh 1 



V4EC 



where a turning point appears for t = 0.) 



It can easily be seen from ( |141| ) to (144) that the harmonic time r is finite for the full 
range — oo < t < oo and goes to constants when \t\ — > +oo. 

For the spatial volume of the universe we have the asymptotic behaviour 



T/~exp j2\C\q\t 



t 



oo, 



(145) 



for all cases (|137|) to (|T4q ). 

In the Euclidean region ( |33|) holds unchanged, since the Wick rotation r — > —it has to 
be accompanied by the transformation a? — > ia 3 (p 3 — » ip 3 ). This means that the parameter 
e in the constraint equation (|29| ) remains unchanged (see 

Thus, the Euclidean metric reads 



ds 2 = dt 2 + a\{t)g^ + ... + a 2 n (t)g( n 



(146) 



where each scale factor has its own turning point at "time" ti, when ^-Oj = 0. The metric 
has its most simple and symmetric form in the case e = {a 1 = 0, % = 1, . . . , n), whence 

ds 2 = dt 2 + (V 1 ) + ...+ 5 ( n )) , (147) 

where V is given by equations ( |137[) to ( |140| ), and the throat is located at t = 0. In the limit 
|t| — > oo, the metric ( 146 ) converges to 



ds 2 = dt 2 + exp 



g (1) + ■■■ + g {n) 



(148) 



which describes an asymptotically anti-de Sitter Euclidean universe. Thus, the metric ( |146| ) 
describes asymptotically anti-de Sitter wormholes. 

The structure of a universe for models with classical Euclidean wormholes is shown 
schematically in Fig. 5 und Fig. 6 for the symmetric case a 1 = (i = l,...,n) with a 



metric (147). (Note that , (|1471 ) has a common scale factor for all factor spaces, while for the 
universe of Eq. (87) in ||23[], sketched by similar looking Fig. 3 resp. Fig. 4 there, internal 
spaces are static, but external space is not.) There are two qualitatively different pictures. 
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The first case (see Fig. 5) takes place for E > (B > 0, C < 0) and for E > (B, C < 0) and 
describes asymptotically an anti-de Sitter wormhole and a baby universe which can branch 
off from this wormhole. The second case (see Fig. 6) takes place for -B 2 /A\C\ < E < 
(B > 0, C < 0) and describes, besides wormhole and baby universe, an additional parent 
instanton which is responsible for the birth of the universe from nothing. 



6. RECONSTRUCTION OF THE POTENTIALS 



The effective perfect fluid Lagrangian ( jl3| ) has its origin in the scalar field Lagrangian (||). 
In this Chapter, we calculate a class of potentials £A )((pW) which ensure the equivalence 
of these Lagrangians. The procedure of potential reconstruction was proposed in and is 
applied in the following. For the integrable 3-component model holds 

( a) _ i r 4^ny)dv (a) 



pW=±.v I v TTv+VO' a = 1,2,3, (149) 

where the energy density p^ is given by (14) and = 2, = 1, = 0. p^ is a 



constant of integration. We should stress that equation ( |149[) was obtained for Lorentzian 
regions. As a result, we get the scalar fields ip^ as a function of the spatial volume V. 
Inverting this expression, we find the spatial volume as a function of the scalar field p^ a \ 
V = V(p^), and consequently, a dependence of the energy density on the scalar field 
tp( a \ p^ = p( a \ip( a ">). Then, using Eqs. (0) to @, we find the potential U^ a \p^) in the 
form 

U {a) {p {a) ) = ~ (2-a (a) ) p {a) (ip {a) ), a =1,2,3, (150) 

where ( ^ 

p (a) =A (a) ly^M)] -0 ° . (151) 



The third component of the scalar field has = 0. Then, from (|150| ) and (151) it follows 
that U^ 3 ', and p^ are constant. This scalar field component with the equation of state 
i-*( 3 ) = — p^ is equivalent to the cosmological constant A = k 2 U^ = k 2 A^> = C. For 
a (1) = 2, we have = (free scalar field). In this case, the scalar field <pW is equivalent 
to a ultra-stiff perfect fluid (P^ = p^). Equation ( |149| ) reads in this case 



<JV - (pW = ± / =, (152 

Y ^° q J VVE + BV + CV 2 V ' 



where E and B are defined by (|S4]) and (^5|) respectively. 
A consequence of (p^) is 



- = ±v / 2 ^(i) r . (153) 

This result is expected for a free minimal coupled scalar field in the harmonic time gauge 
where <p^ = 0. After integration in ( |152D , 



(i) (i) ■ g VBV + CV 2 
^>-<Po =Ti— > E = > ( 154 ) 



{1) (1 ) 2E + BV-2VER 
pU_ v U=±__] 11 _ , £>0, (155) 
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I AW i 



BV-2\E\ 



(1) = ±i-^=^arcsin^^4^, E < 0, B 2 - 4EC > 0, (156) 



with fi:=£ + W+ CV 2 and |A| = B 2 - 4EC. For E < (i.e. < 0) this scalar field 
is imaginary. 

Let us now consider the second component with = 1. The scalar field is equiva- 
lent to dust (P (2) = 0). Equation ( p49|) reads now 



(157) 



9?( 2 ) is imaginary for < 0, i.e. B < 0. The integral in (|157| ) is an elliptic one and, in 
general, it is not possible to express it by elementary functions. But in the particular case 
E = 0, which expresses the asymptotic behaviour of the scalar field (|157|), we get 



(2)_ M ( 2 ) =T ^ ai( , )(h ( | + _ r 



n 



V2 



<p&) - ip {2) = =p— artanh ( 1 
Kq 

(2) (2) -^-V 72 " , 

ip y J — (fr. = =fi — arctan 
Kq 



C_ 

\B\ 



C \V3 
B 

1/2 
1/2 



B 



V 



V-l 



the volume of the universe 



V = — sinh" 
C 

y = ^cosh- 

K = - 1 — 1 cos ' 
c 



^ (V 2 ) 
^ r„( 2 ) 



(2) 



V2 



K 1 ■ ( (2) (2) 



v 7 ^ 



and the potential of the scalar field 





c 


sinh 2 


" Kq 
.72 




\C\ 
2k 2 


cosh 2 


' Kq 

.71 




C 

~ 2^2 


cos 2 


Kq 
y/2 



(2) 

(2) 
¥>0 



£,C > 0, 

B > 0, C < 0, 

S < 0, C > 0, 

B,C > 0, 

5 > 0, C < 0, 

5 < 0, C > 0, 

5,C > 0, 

£ > 0, C < 0, 

B < 0, C > 0. 



(158) 
(159) 
(160) 

(161) 
(162) 
(163) 

(164) 
(165) 
(166) 



It follows from these equations that for B,C > and B,C < the volume goes to infinity 
like 



V 



1 



(2)1 



+ 00, 



(2) 



0. 



(167) 



The general expression (157) should have the same asymptotic behaviour in all the cases 
where the limit V — > +oo is permitted, because we can drop in this limit the term E in the 
denominator of (|157j). 



If E > 0, from (fL57|) results 



(2) (2) 



2 AM /E 



(168) 
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Let us now consider two particular cases of ( |157| ) for E ^ 0. The first case is that one when 

(1 2) 

the classical trajectory has two turning points Vq , i. e. when either B > 0, C < or 
B < 0, C > ( for both the cases B 2 > AEC). Then (see equation (3.131) in [JO]), 



p(2) = ±A ^ ^F(V>|m) (169) 
K( l J\B\ + J\K\ 



where F{ip\m) is the elliptic integral of the first kind [42| and 



4> = arcsin ^ (F {2) - V)/{V^ 2) - V {1) ), F (2) > V > F (1) , 5 > 0, C < 0, (170) 
if) = arcsin yJv/V^\ V^ 2) > >V, B < 0, C > 0, (171) 



z/> = arcsin yJ(V -V$ 2) )/(V -V^), V > V^ 2) > B < 0, C > 0, (172) 



m 



> V > Vr 



(1) 



! 



B > 0, C < 0, 



m 



F (1) /^ (2) , V {2) > F (1) > V, B < 0, C> 0, 



v > y (2) > f (1) , b < o, c > o. 



The turning points are 



(1,2) 



\B\ T# 



2|C| 



(173) 
(174) 
(175) 

(176) 



The minus sign is related to Vq , the plus sign to V$ ', and |A| = B — AEC. The scalar 
field (p( 2 ^ is imaginary for B < 0. 



(2) 



With the Jacobian elliptic functions [42], inverting ( |16S| ), the volume of the universe is 
given by 



V=Vr 



(2) 



sn^ 



(2) 



V = sn z V c 



vf } ), 

2t/(1) 
> 



v£ 2) >V> V (i), B > 0, C < 0, 
F (2) > F (1) > V, B < 0, C > 0, 



(2) 



2t/(1) 



sn^K 



1 — sn 2 



v > vi 2) > \ , 



! 



B < 0, C > 0, 



(177) 
(178) 

(179) 



where sn = sn {l\ip^ \mj = sinip and I 1 1 



2^ 



The corresponding potential 

Kqy/\B\+\A\ 

terms are then given as 

jj{2) _ ^(2)/2V (see (pa and ( pH|) for = 1). According to the 
properties of the Jacobian elliptic functions |42[ [ asymptotic estimates for ( 178 ) and ( |179| ) 
are V « (q 2 E/2AW)(ipW) 2 for (t^ 2 )] -> (in accordance with (MB)) and V ~ l/l^ 2 ) [Mbr 



(2)1 



(in accordance with (|l67p ) 



Another particular case for C > is that with E > B 2 /4C. Here (see (3.138 (7)) in pjl), 
we obtain 

1 



<p {2) = ±*^jE]b^ mm) ' (180) 



where 



and 



nq 



ip = 2arctan dV/JE/C 



m = J {2VEC - B) /4VEC. 



(181) 
(182) 
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Inverting equation (180), we obtain 



sn z v 



sn^ 



sn^ 



E/C] 



E 
C 



(183) 



with sn 



sn(l2(f^\m) = simp and I 2 
1^(2) 1 2 — * 00 f° r \ ~^ (i n accordance with ( |l6 

sign, V « (g 2 £;/2^( 2 ))(v?( 2 )) 2 -> for |</?( 2 )| -> (in accordance with ( g68| 
scalar field potential, we have to substitute (|l83|) into U&\ipW) = A^/2V. 



V 



i 



± J^EIQ^ rji ■ For the branch with the plus sign, 



and for the branch with the minus 
To find the 



7. SOLUTIONS TO THE QUANTIZED MODEL 



At the quantum level, the constraint equation fl25| ) is replaced by the Wheeler-DeWitt 
(WDW) equation. The WDW equation is covariant with respect to gauge as well as minisu- 
perspace coordinate transformations [25|. In the harmonic time gauge [13, |25l] it reads 



i d 2 d 2 



2dz ' 



^E^ + ^E A(a) «p(fc< tt V) ] * = 0. 



(184) 



0=1 



Formally, this WDW equation has the same structure as that of |23| |. However, on the 
semiclassical level the dynamics of the universe is quite different for the models in both the 
papers. Semiclassical equations were considered in [43]. 

We look for solutions of ( |184[) by separation of the variables and try the ansatz 



■$r(z) = $(z°)exp (ip • z) , 



(185) 



where p := (p , . . . ,p n ) is a constant vector, z := (z , . . . , z n ), Pi = p % and p • z : 
TlClvd- Substitution of fll85p into (p^ ) gives 

72 -i n-l 

2dz° 2 ' 2 



1 d 2 



-, n—1 m 

+ -^2 (pi? + ^ 2 E A(a) ex p ( k(a) Q z ° 



i=l a=l 

For the integrable 3-component model this equation reduces to 

1 d 2 



$ = 0. 



(186) 



$ = E$ 



(187) 



in the notation of (gj) and (p5|). Following [23], we rewrite this equation like 

d 2 $ d<S> 



dx 2 dx 



E + Bx + Cx z 



$ = 0, 



(188) 



where E = 2E/q 2 , B = 2B/q 2 , C = 2C/q 2 , and x = exp{qz°) (x is identical to the spatial 
volume of the universe V) . (fL88|) is equivalent to the Whittaker equation 



di 2 



+ 



1 B/T £ + 1/4 
"4 + -^- + — ^2— 



y = o, 



(189) 
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where T := ±2iV~C, £ := Tx, and <I> =: x 1//2 y(£) and also equivalent to the Kummer 
equation 



d w , ^dw 



1 



B 



2 +/i T 



w = 0, 



(190) 



where fi 2 := —E and =: x l l 2 exp f — |c) £ 2+/ ^(£)- ^ n the f* rs t case, the solutions are the 

Wittaker functions @ yi := M fc]jU (£) and y 2 := Wfe )( u(£) with k:= B/T and /i 2 := In the 
second case, the solutions are the Kummer functions iui := M(a, b, £) and u>2 := 17 (a, 6, £) 
with a := i + /u — B/T and 6 := 1 + 2/^. 

The general solution of equation ( |184| ) for the 3-component model is 



(191) 



= E / d n ~ 1 pa( P )exp(fp-z) $f (exp(^ )), 

„■ 1 O ** 



i=l,2 



where < I ) ^' 2 ' 1 = 77=yi,2(£)i or ^e'^ = 77J ex P ( — |£j £ 2+Mu, i,2(£)- It is convenient to set 

T = +2iVd for C > and T = -2iVC for C < 0, and p := + v / ^E. 

In [43], it was argued that the parameter E can be interpreted as energy. So, the state 
with E = 0, vanishing momenta pi [i = 1, . . . ,n — 1), and = (absence of free scalar 
field excitations) is the ground state of the system. Thus, its wave function reads 



(exp(gz°) 



1,2. 



(192) 



The limit of large spatial geometries in (|188|) z° — » +oo (remember x = V = exp(qz )) is 
equivalent to B —> 0. In this limit, the Wittaker functions reduce to Bessel functions |42]| , 
namely 



for C > and 



k-*0 



VI,, 



VK„ 



\C\V) , 
^\V 



(193) 
(194) 

(195) 
(196) 



for C < 0. 



Following the ideas of [43, 44], one can demonstrate that for C > the wave function 

' /2C ' 



fe— >o 



-V 



~ cos 1/ 

V^oo \ q 



corresponds to the Hartle-Hawking boundary condition [45] and the wave function 



v 



(2) 



H, 



(2) 



/ 2C, 



-V 



~ exp 

V-+oc 



(197) 



(198) 



corresponds to the Vilenkin boundary condition 
In the case C < 0, we get 



X J,HH _ k r 



V2C 



V" 



~ exp 



\J2C 



V 



(199) 
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The potential has for B > 0, C < a well and for E < the energy spectrum is discrete 
(see Fig. 2). In this case, the finite solutions of the wave equation ( |187| ) | |47| are 



<D n = exp(--£)^M(-n,&,£). 



(200) 



The energy levels are given by 

-E, 



B 



2y/\C\ V2 

n is a non negative integer and restricted to 



n + 



n < 



B 



qVW\ 



(201) 



(202) 



Thus, the discrete spectrum has a finite number of eigenvalues. If — §= < i, there is no 

discrete spectrum. It was demonstrated in [23] that the wave functions ( |200D satisfy the 
quantum wormhole boundary conditions Eq]. 



8. CONCLUSIONS 

We considered the generalization of a homogeneous cosmological model of Bianchi type I to 
an anisotropic multidimensional one with n > 2 Ricci-flat spaces of arbitrary dimensions, 
in the presence of m homogeneous non-interacting minimally coupled scalar fields. Under 
certain conditions these models are equivalent to multidimensional cosmological models in 
the presence of an m-component perfect fluid with equations of state = — l^j 

with matter constants for a = l,...,m. Using this equivalence, for m = 3, we find 
integrable models when one of the scalar fields is equivalent to an ultra-stiff perfect fluid 
component, the second one corresponds to dust, and the third one is equivalent to a vacuum 
component. Dynamics of the universe was investigated in general, as well as in a particular 3- 
component integrable case. For integrable models, there are four qualitatively different types 
of evolution of the universe, depending on the potential U(zq) (see Fig. 1 and Fig. 2), but in 
all four cases the universe has a Kasner-like behaviour near the cosmological singularity (see 
|123| ). In the cases where the universe can expand to infinity, an isotropization takes place 
which results in an asymptotically de Sitter universe (see ( [L33D and (1134 )). 

In quantum cosmology, instantons, solutions of the classical Einstein equations in Eu- 
clidean space, play an important role, giving significant contributions to the path integral. 
They are connected with the changing geometry of the model. We found here three inter- 
esting types of instantons. The first one (see ( |107| ) and ( |109| )) describes tunnelling between 
a Kasner-like universe and an asymptotically de Sitter universe. Sewing a number of these 
instantons (see Fig. 4) may provide the Coleman mechanism for the vanishing of the cosmo- 
logical constant. Another type of instanton (see (|82]) and ( |108| )) is responsible for the birth 
of the universe from " nothing" . It was shown that corresponding Lorentzian solutions (|8l| ) 
and (|105| )) can ensure inflation of the external space (see ( p.36|) ) and compactification of the 
internal ones. This problem needs a more detailed investigation in a separate paper. The 
third type of instantons (see equations ( |137| ) to ( 140[ )) describes the Euclidean space which 



has an asymptotically anti-de Sitter wormhole geometry (see Fig. 5 and Fig. 6). 
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The scalar field potentials U {a) (^ a) ) {a = 1, . . . , m) can be reconstructed by the method 
described in We performed this procedure for integrable models, and exact forms of 

potentials were presented in section 6. 

The equivalence between a scalar field and a perfect fluid component helps also to inves- 
tigate the quantum behaviour of the universe. We obtained the Wheeler-de Witt equation 
from the effective perfect fluid Lagrangian. Exact solutions are found, some of which describe 
cosmological transitions with a signature change of the metric, e.g. universe nucleation as 
quantum tunnelling from an Euclidean region. Other solutions are given as quantum worm- 
holes with discrete spectrum (see ( |201| )). 
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FIGURE CAPTIONS 



Fig. 1. The potential Uq(z°) (solid line) and the energy levels E (dashed lines) in the cases 
B, C > and B, C < 0. Lorentzian regions exist for E > Uq(z°). 



Fig. 2. The potential Uq(z ) (solid line) and the energy levels E (dashed lines) in the cases 
B > 0, C < and B < 0, C > 0. In the former case we get a potential well, and in the 
latter case we obtain a potential barrier. Lorentzian regions exist for E > Uq(z°). Here, 
U m = B 2 /AC and z° m = \ In \B/2C\. 



Fig. 3. The qualitative shape of the instanton ( |107|) (or ( |109|) ). The instanton describes 
tunnelling between a Kasner-like (baby) universe and a de Sitter universe. 



Fig. 4. Examples of instantons constructed by sewing together several copies of the the 
instanton illustrated in Fig. 3. Such instantons may describe tunnelling between (a) Kasner 
and de Sitter universe, (b) two Kasner universes, (c) two de Sitter universes. 



Fig. 5. An asymptotically anti-de Sitter wormhole is shown schematically for energies E > 
(B > 0, C < 0) and E > (B, C < 0) in the symmetrical case a* = (i = 1, . . . , n). (Note 



that due to suppressed internal degrees of freedom this looks like Fig. 3 in ^3| , although the 
solution and its context are quite different there.) 



Fig. 6. The qualitative structure of the universe is is shown schematically for energies 
— U m < E < (B > 0, C < 0) in the symmetrical case a 1 = (i = 1, . . . , n), (Note that due 
to suppressed internal degrees of freedom this looks like Fig. 4 in |P3| 1, although the solution 
and its context are quite different there.) 
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